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We study the classical anisotropic ferromagnetic spin chain with frustration. The behavior of 
soliton and kink solutions in the vicinity of the ground state phase transition from the ferromagnetic 
to the spiral phase is studied. The dependence of the soliton energy on small anisotropy parameter is 
established using scaling estimates and numerical minimization of the energy functional. Conditions 
of the existence of the solitons are determined. It is shown that solitons survive in the spiral phase 
, though with some restrictions on their size. A comparison of the energies of the classical solitons 

and the bound magnon complexes in the quantum model shows the functional similarity between 
them. The influence of the finite-size effects on the soliton states is studied and it is shown that the 
£SJ \ localized solitons originate from the uniform state when the system size exceeds some critical value 

depending on the anisotropy. 

< 

I. INTRODUCTION 

Lately, there has been considerable interest in low-dimensional spin systems that exhibit frustration A very 
interesting class of such systems is chain compounds consisting of edge-sharing CuO± units. Recently, a variety of 
these copper oxides were synthesized and found to show unique physical properties @, [H, 0, 0] ■ The frustration 
in these compounds arises from the competition of exchange interactions between magnetic Cu 2+ ions carrying spins 
& ' 1/2. Due to a specific geometry of these systems (Cu — O — Cu angle is close to 90°) the nearest-neighbor (NN) 
interaction is ferromagnetic, while the next-nearest-neighbor (NNN) one is antiferromagnetic and absolute values of 
these interactions are comparable 0, 0j- An appropriate model describing the magnetic properties of such copper 
oxides is so-called F-AF spin chain model the Hamiltonian of which has a form 

c : 

O . N N 

L^£j| H = J\ y^ y (^n^n+l + Sn^n+l + ^l^n^n+l) + ^2 (fin^n+a + + A 2 S*S* + 1 ) (1) 

n—l n—l 

t-H ■ 

J> ■ where J\ < and J2 > 0. 

This model is characterized by a frustration parameter A = Ji/\Ji\- The quantum F-AF s = 1/2 model has been 
intensively studied last years [1, [t| EH EH, E2> E!) EH EE EH • Most of studies of this model are related to the isotropic 
case (Ai = A 2 = 1). It is known that the ground state of the isotropic model is ferromagnetic for A < 1/4. At A = 1/4 
, the phase transition to the incommensurate singlet phase with spiral spin correlations takes place. Remarkably, this 
transition point does not depend on spin value s at N — > 00. It was shown also that the F-AF s = 1/2 model 
with anisotropic interactions has a rich phase diagram 17]. In our papers [HI. fli^| we investigated weakly anisotropic 
(Ai > 1, A2 = 1) quantum spin-1/2 model fTJ). H was shown that even small anisotropy essentially effects on the 
properties of the model. In particular, the transition point from the ferromagnetic to the spiral- like ground state 

>! shifts from A = 1/4. 

V^J ' An interesting feature of the quantum anisotropic model ([1]) is the existence of the multimagnon bound complexes 
in the ferromagnetic phase. These complexes govern the low-temperature thermodynamics (l9l.l20j. It was noted also 
that the NNN interaction strongly affects the excitation spectrum especially in the vicinity of the transition point 
A = 1/4. 

It is known that there is a close relation between the multimagnon bound complexes in the quantum spin models 
and soliton and kink excitations in the classical counterparts. In particular, these soliton states have been studied 
extensively for the classical easy-axis ferromagnetic chain and a connection between them and the bound magnon 
complexes in the quantum s = 1/2 model ([1]) at J 2 = was discussed [21] . On the other hand, the influence of the 
frustration on the solitons has not been considered before. The F-AF model represents a suitable model to study this 
problem and to investigate the connection between the quantum excitation spectrum and the soliton solutions of the 
classical frustrated spin model. 
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In [18[ we showed that the ground state phase diagram of the quantum model (|TJ) with both Ai =^ 1 and A 2 ^ 1 are 
qualitatively similar to that for the model with the anisotropy of the NN interaction only. Therefore, for simplicity 
we consider model (p} with Ai = A > 1 and A2 = 1. In this case model (JT]) takes the form 

H = - (S:S* +1 + SyS v n+1 + - A, 2 ) + A^(S„ • S„ +2 - s 2 ) (2) 

where we put | J\ | as an energy unit and added constant shifts to secure the energy of the ferromagnetic state to be 
zero. 

The phase diagram of the classical F-AF model with A > 1 consists of the ferromagnetic and the spiral phases. 
The ferromagnetic ground state is simple, i.e., it has all spins parallel to the Z axis. However, the soliton excitations 
in this phase are not trivial especially near the transition point between the phases. Our main goal is to study the 
behavior of the solitons in the vicinity of the isotropic transition (IT) point (A = 1, A = 1/4) and to compare it with 
the properties of the excitations of the quantum model. As it will be shown the frustration effects strongly modify 
the soliton states especially for small anisotropy. In particular, the exponent characterizing the power dependence of 
the gap (the soliton energy with respect to the ground state) on the anisotropy is different from that for A = 0. 

The paper is organized as follows. In Sec. II we represent the known results for soliton solutions of the anisotropic 
ferromagnetic chain (model @ at A = 0). In Sec. Ill we consider the classical continuum F-AF model in the vicinity 
of the ground state phase transition from the ferromagnetic to the spiral phase. We deduce the corresponding energy 
functional and establish the scaling form of the soliton energy as a function of the anisotropy and the frustration 
parameter. In this section we also obtain asymptotes of the soliton solutions at large distances and determine the 
necessary conditions of the soliton stability. In Sec. IV we present results of the numerical minimization of the energy 
functional which confirm the scaling estimations. In Sec.V we study the finite-size effects on the soliton solution. In 
Sec. VI we show that the behavior of both classic solitons and m-magnon quantum spin excitations are functionally 
similar to m-boson bound complexes of the Bose model with the attractive interaction if m is not large. In Sec. VII 
we give a summary of results. 



II. CLASSICAL MODEL FOR A = CASE 



In the classical approximation the spin operators S n are replaced by the classical vectors of the fixed length s 
which are parameterized by spherical coordinates 

~£) n = s(cos Lp n sin0„, sin ip n sm9 n , cos 8 n ) (3) 
In terms of the angles 9 n and ip n the discrete classical F-AF model <j2j) takes the form 

E = fi2 U 1 ~ cos ( e n+i - °n)] - A[l - cos(0„ +2 - 6 n )] + a[l - cos(0 n+1 ) cos(0„)]} 



+s- 



y^{[l - co$(cp n+1 - ip n )] sin6>„ sm8 n+1 - A[l - cos(<^„ +2 - <£„)] sin6>„ sin6>„ +2 } 



(4) 



where a = A — 1 > 0. 

In this Section we briefly review the known results for the classical Heisenberg ferromagnetic chain with an easy-axis 
anisotropy, i.e., model (HJ with A = 0. In this case the model is exactly solved in the continuum limit [22I [23j] . 

The continuum limit of the classical model assumes that the vectors can be replaced by the classical vector 
field ~£>(x,t) with slowly varying orientations, so that 



^ n+1 O r 



dx 



(5) 



where the lattice constant is chosen as unit length. The direction of the vector field 1>(x,t) is determined by two 
angular variables <p(x,t) and 8(x,t) according to Eq.([3]). 

The dynamics of the vector field ~3(x,t) is governed by the Landau-Lifshitz equation: 

a * " - (e) 



Of 



SE 



where we put % = 1. Here E is the energy as a functional of the vector field S (x,t). Using the continuum approxi- 
mation |5]), Hamiltonian ([2|) goes over into the well-known energy functional: 



E 



dx 



1 (djf 

2 I ~dx~ 



+ a (, 2 - Si) 



(J) 
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Classical equation of motion ([6]) for model ([7]) has two constants of the motion: the magnetization (continuum 
analog to a number of magnons) 



M 



tj da;(l-cos0) (8) 



where we subtract a constant to make M finite for solitons, and the momentum 

P = s J dx(l -cos0)|| (9) 

The ground state of ([7]) is the ferromagnetic configuration with all spins parallel (or antiparallel) to the Z axis, i.e. 
9 = (9 = 7r). The soliton configurations are specified by the boundary conditions 9 — > (9 — > tt) at x — > ±oo. 

Fortunately, the classical equations of motion (jSJ) for model (J7J) are exactly solvable (24|. In particular, the energy 
for a given values of M and P is 



., — cosh(MV2a/2s) - cos(P/2s) 



E M p(a) = isW2a - — '—}= — — - (10) 

1 sinh(MV2a/2s) V ' 

Remarkably, Eq. (flQ|) reproduces the exact result for the energy at a <C 1 of the M-magnon bound state [25| for the 
most quantum case s = 1/2. 

Here we will show that simple scaling arguments are able to establish the correct scaling dependence for the energy. 
For this aim we rescale the coordinate x = ^/a 1 ^ 2 and introduce the normalized spin vector field ~rt(x, t) = S(x, t) / s 
in Eq.([7]), which results in 



E = s 2 a 1 ' 2 J d£ 



1 / &rt x 



nl) 



(11) 



We notice that the integrand in equation (|11|) does not depend on any parameters and is expressed through the 
normalized spin vector field. This means that the energy scales as E ~ s 2 a 1 ^ 2 , which is correct for a kink or for a 
large soliton excitations (see Eq. (|10|) ). However, we can make one more step. The same procedure for magnetization 
© gives 

M = sa~ 1/2 [ d£(l-n z ) (12) 



This expression means that the magnetization forms a scaling parameter Mo 1 / 2 / s. In a similar manner one can 
find that momentum (|9|) produces a dimensionless parameter P/s independent of a. 

Thus, the energy of a soliton of the size M having a momentum P can be written in a form: 

E = s 2 a 1/2 f{Ma 1 / 2 /s, P/s) (13) 

where / is a scaling function, which can not be found in the framework of this scaling estimate. Fortunately, for model 
this function is known exactly (fTT)|) . which validates the above scaling arguments. Thus, simple scaling estimates 
allowed us to establish the scaling parameters and the scaling form for the soliton energy. 

There is one more important fact which is worth noting here. The exact static solution for a kink in the discrete 
XXZ model with classical spins has been constructed by Gochev in Ref.[26[. He found that the kink energy has a 
form 

£ kink = 2sVA 2 - 1 (14) 

which reproduces the exact result for the case s = 1/2 [25| . This allows us to assume that the energy of kink (or large 
soliton, which has double energy of kink) for XXZ chain with any value of s has an universal dependence on A ([Tl 
We will see that this universality is destroyed in case A ^ 0. 



III. CLASSICAL CONTINUUM SPIN MODEL NEAR THE IT POINT 



The NNN term in model ([2|) causes the frustration in the system and immediately destroys the integrability of the 
model. However, the behavior of the system for < A < 1/4 and weak easy-axis anisotropy remains very similar 
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to the case A = 0. In the classical and continuum approximation the Hamiltonian reduces to an energy functional 
([7]) with renormalized factor (1 — 4A)/2 at (dS / dx) 2 . In particular, this means that the found scaling form for the 
soliton energy (| 13j) remains valid. 

The situation drastically changes near the IT point (A = 1/4, A = 1) where the ground state phase transition in the 
isotropic case takes place. Let us deduce the energy functional describing the vicinity of this point. The continuum 
approximation assumes that the differences (0 n +i — On) an d {<Pn+l ~ <Pn) are small and Eq. (J4j) can be expanded in 
these differences. However, in the vicinity of the IT point it is more clear and instructive to derive the continuum 
approach in terms of a spin vector field. For this aim we rewrite Hamiltonian ^ in the form: 

H = \ £(S n+1 - 2S„ + S n _!) 2 - ! ]T(S„ +1 - S^O 2 + a^> 2 - S z n S z n+1 ) (15) 

where the parameters a = A — 1 and 7 = 4 A — 1 are small in the vicinity of the IT point. 

In the classical and the continuum approach the expression containing spin operators in the first term in Eq. (|15p is 
replaced by 

1? q-<? _l"c? ~ d 2 !)(x n ) 

<J n+l ~ ^ <~> ri + O n -l 7^ 

and that in the second term in Eq. (fl"5)) according to Eq. ([5]) . 
Thus, Hamiltonian (|15[) is mapped to the energy functional 



E = / da 



1 (d 2 T 
8 I dx 2 




S 2 Z ) 



(17) 



The effect of the fourth-order term (d 2 ~§ /dx 2 ) 2 in the energy functional like (fT7|) was studied before [13], but it 
was considered as a small correction to the main contribution given by the term (d~§ / dx) 2 . On the contrary for our 
model (fl~7|) near the IT point the fourth-order term (d 2 ~§ / dx 2 ) 2 becomes the leading one. The appearance of the 
fourth-order term is related to the fact that the one-magnon spectrum in the IT point becomes e(k) ~ k 4 . 

The energy functional (|17p in terms of the angular variables <p(x) and 9{x) has a rather cumbersome form: 



1 sm 2 I 



E = S - J dx [6" 2 + 9' 4 ~ 4 7 #' 2 + y 2 + v 14 - 4 7 (^ 2 + 8a) 1 

+ ^ J dx [(4 - 2 sin 2 9)<p' 2 6' 2 + {2ip"<p'6' - Lp' 2 9") sin(26»)] (18) 

where the prime denotes the space derivatives d/dx. One can check that Eq.lfTS]) represents the leading terms in the 
expansion of Eq.(f4]) in small differences (9 n +i ~ 9 n ) and (f n +i — fn)j which is actually assumed in the continuum 
approximation. 

At first, let us study the phase diagram of the classical continuum model (fl~8|) . Numerical calculations (details will 
be given below) confirmed a natural assumption that the ground state spin configuration for the easy-axis anisotropy 
case (a > 0) corresponds to the choice <p = const. In this case the energy functional (JTSJ) simplifies to 



E = — J dx(9' u + 6»' 4 - 4 7 <r + 8a sin" 9) (19) 

Variation of the energy functional (fT!)|) in 9(x) leads to the Euler equation 

\q»» _ Iffttf* + ld " + a sm{29) = (20) 

The ground state in the ferromagnetic phase has a trivial solution 9 = (or 9 = n) with zero energy. In the 
isotropic case (a = 0) the transition from the ferromagnetic to the spiral phase takes place at A = 1/4. For A > 1/4 
the solution 9(x) has a pure spiral form 

9 sp (x) = ±^x (21) 
which evidently has zero mean spin projection on any axis (S x VtZ ) = and describes the spiral in the XZ plane. 
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FIG. 1: The phase diagram of the classical F-AF model. Thick line is the phase boundary between the ferromagnetic and the 
spiral phases. Dashed lines are boundaries between regions of different asymptotic regimes for large static solitons. 



In the anisotropic case the spiral function (|2Tj) does not satisfy Eq. (j20]) . However, it represents a good approximation 
for the solution in the spiral phase. Using this function as a variational one for the energy functional (|19p we find 
that the transition between the ferromagnetic and the spiral phases occurs on the line 7 = ^fu. 

Assuming that the correction to the spiral solution 9 sp (x) is small, we found from Eq. (f2"0|) that the first correction 
has the oscillating form: 

6{x) « y^x - ^ sin^x) (22) 

This function is not a pure spiral: due to the oscillating correction the spin vectors prefer to direct along the Z 
axis, so that (S^) > (S 2 + S^). However the mean values of all total projections remain zero: (S x .y : z) — 0- 
The calculation of the energy functional (fT9|) with the function ([22]) yields 

E s P _ 7 2 - a Q 2 

s 2 N ~ 2 128 7 2 1 ' 

One can see that the correction to the energy is really small in the spiral phase region. Energy (f2"51) becomes negative 
for 7 > 0.992 v / a. This means that taking the correction in Eq. ([22|) into account shifts the transition line from the 
ferromagnetic to the spiral phase on the value less than 1% (the smallness of the correction in Eq. (|22"|) allows us to 
leave the name spiral for this phase). The phase transition from the ferromagnetic to the spiral phase at 7 w 0.992 v / a 
is obviously the first order one. The transition line 7 = Q.QQZ^Ja is shown in Fig[TJ 

It is worth noting here that the found behavior of the transition line corresponds to the classical version of spin 
model @ and it is shifted for the quantum spin-s model. Quantum corrections for s > 1 change jnumerical coefficient 
in Eq. (|2"5|) . But for the s = 1/2 case they lead to another form of the transition line: 7 ~ cv 3 / 5 .18]. 

Now let us study the excitations of model (TP?]) . Similar to the case A = the lowest configuration in the ferro- 
magnetic region for a given (S z ) is described by the static soliton-like solutions of Eq. (|20|) with ip = const over the 
ferromagnetic configuration 9 = (or 9 = ir). Unfortunately, we could not find the exact soliton solution near the IT 
point. However, it is possible to determine the scaling dependence of the energy of soliton of size M on parameters 
a and 7 using scaling estimates as we did for the case A = 0. Near the IT point we perform the rescaling x = ^/cv 1 / 4 
and S(x,t) = slt(x,t) which transforms the energy functional (fT"7|) to 



E = ,s 2 a 3 / 4 J df 



1 / d 2 Tt\ [i ( dlt 



8 V <9£ 2 / 2 V <9£ 



fl - nl) 



(24) 



with fj, = j/*/a. 

Since the integrand in Eq. depends on the parameter /i only, we conclude that the energy of a kink or of a large 
soliton is E = s 2 a 3 / 4 f(/j). As a result of rescaling the magnetization (jSJ) becomes: 



M = sa- 1 ^ 4 I d£(l - n z ) (25) 
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So, the magnetization produces a scaling parameter v = a^l^Mjs. The momentum © forms the same parameter 
k = P/s as for the case A = 0. 

Thus, the energy of a soliton of the size M can be represented in a form: 

E = s 2 a 3 / 4 f(»,v,k) (26) 

where the scaling function / can be found numerically. This scaling form is exactly the same as was found for the 
quantum spin- 1/2 model near the IT point (l9j |. 

We note that similar scaling procedure for the D-dimensional version of model (|17p gives the scaling form for 
the energy of the soliton of the size M (even small-amplitude solitons are stable for model (fTT)) in two- and three- 
dimensional cases [28[): 

E = s 2 a 1/2 f(^,a 1/2 M/s,k), D = 2 

E = s 2 a l ' i f{n,o?' i M/s,k), D = 3 (27) 

and the case D = 4 requires special treatment. 

Though the exact solution of the corresponding equation of motion ([6]) is unknown it is possible to identify the 
necessary conditions for a stability of solitons and establish their asymptotic behavior. It is convenient to introduce 
the complex function 

ip = n x + in y (28) 

so that ip — e llp sin# and n 2 — 1 — \ip\ 2 . 

Equation of motion ([6]) for -0 -r) near the IT point reads 

dip 1 d A ip 1 , d A n z d 2 ip , d 2 n z 

i ^ = r*w _ t^w + Mnz w " ^ «f + ^ (29) 

where r = sa 3 ^ 4 t is rescaled time. 

Solitons are localized objects and, therefore, it requires that j^l — > far from the soliton center. So, at large 
distance from the center one can linearize Eq. (|29p in ip by putting n z = 1: 



.dijj _ 1 <9V d 2 ijj 



i-5Z = 7«^+^ + 2tf (30) 



We seek the asymptotic of ^> in conventional exponential form 

T ) = ex P [—iur — k(£ — vt)] (31) 

where v and u> are normalized linear and angular velocities defined as 

1 dE _ df{n, v, k) 
as dM dv 
_ 1 df(n,v,k) 

~ aV^s dk (62) 
Substituting Eq. ([3T|l into Eq. ([3l]|) we obtain equation for k 

K 4 

— + [in 2 - ivn + 2 - uj = (33) 

This equation has four roots. The condition of the soliton existence is that all roots have non-zero real part, 
which secures the decay of solutions at |£ — vt\ — > oo. The region in (v,uj) plane where ip exponentially vanishes at 
|^ — vt\ — > oo is defined as lj < loq(v). The dependence u>o(v) is parametrically definable function 

v = q 3 ~ 2fiq 

uj q = 2 + M9 2 -|9 4 (34) 
where the parameter q runs from — oo to oo for /i < and \q\ > for /i > 0. 
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FIG. 2: Boundaries of the existence regions for soliton in (v,co) plane for different fi. Allowable values of v and to lie below the 
boundaries. 

Equations (|34| of the boundary of the soliton existence region coincide with the condition of the spin- wave instability. 
The existence regions for a few values of the parameter fi are shown in FigJ^l For fi < the region has a form similar 
to the case A = with the quadratic dependence uiq(v) near the maximum point (v = 0, to = 2). For the case fi = 
this dependence becomes loq = 2 — 3t> 4 / 3 /4. When /i > 0, the dependence u>q(v) shifts down and a cusp appears at 
the maximal point (v = 0, to = 2 — /i 2 ). 

The region of the soliton existence contains allowable values of v and ui for a given /i. As it can be seen in Figf5]the 
soliton can exist in the spiral phase (fj, > 0.992). But, for example, the kink is unstable for \i > y2 since the point 
(v = 0, lu = 0) corresponding to the kink does not belong to the existence region. It is worth noting that one should 
be careful in treating of the part of the existence region with large negative lo, where the derivatives of 6(x) and p(x) 
near the center of soliton become large and the continuum approximations Eqs.([5]l. (j!6p become doubtful. 

Within the existence region the roots of Eq. (J33J) have the forms 

ki,3 = ±oi + ib 

K2,4 = ±02 — ib (35) 

where ax, a<i and b are real nonnegative quantities and a\ < a-i- Certainly, one should take into account only the roots 
providing the decay of \i\>\ at the corresponding limit. For determinacy we will consider the limit (£ — vt) — » oo and, 
therefore, analyze the roots Kx = ax + ib and K2 — ci2 — ib. At infinity only the root Kx with the smaller real part is 
vital. So, the asymptotes of the angles at (£ — vt) — + oo are 9 ~ exp[— ai(^ — vt)] and ip = —lot — 6(f — vt). 

The case of the static soliton (v = 0) requires special treatment, because in this case both roots Kx and k-i have 
equal real parts a\ = ai (they immediately split if v ^ 0). Fortunately, for the static soliton the solution of Eq. (|3"3")l 
can be obtained explicitly and this allows to analyze in detail the behavior of the soliton asymptotes for given values 
of n and to. The decaying roots at (£ — vt) — > oo are 

Kx = - 2 \JlP- -2 + LO (36) 

and K2 = k\. 

Simple analysis of Eq. (f3"6"]l shows that there are three regimes of the asymptotic (|3"Tj) . The exponent Kx is real for 
fi < —\/2 — to, which means that 9(£) smoothly tends to zero at infinity. This region lies totally in the ferromagnetic 
phase. In the region —a/2 — lo < /i < \J2 — to the parameter k\ contains both real and imaginary parts. Here 
the superposition of asymptotes (f3"Tj) with Ki and k 2 — k,\ results in the decay with oscillation at infinity: (9(£) ~ 
e -aiC sm (5£) an d ip = const. Finally, for [i > \J1 — lo the parameter k\ is pure imaginary, which implies that the 
solution oscillates but does not decay at infinity. Therefore, in this region there are no static soliton-like excitations. 

As was noted above only one asymptotic with K\ having smaller real part survives at infinity for u^0. This implies 
that the asymptotic of the function 9(£) has no zeros. However, in the region close to the soliton center the function 
can have a finite number of zeros and it is really observed numerically. 

Generally, the relation between the parameters to, v and v, k is unknown analytically. However, it can be found for 
large values of v. For large solitons (v ^> 1) the soliton energy and, therefore, the scaling function /(//, v, k) saturates 
to some finite value. This means that for v ^S> 1 both lo — > and v — > 0. The boundaries between the regions 
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corresponding to different asymptotic regimes for large solitons are obtained from the above analysis and they are 
shown in Fig[TJ As one can see in Fig[T]the oscillating large solitons exist not only in the ferromagnetic phase but, 
partly, in the spiral phase region 0.992 < \x < y/2. The amplitude of the soliton oscillations grow as \i increases and 
the oscillations are most strong in the spiral phase. However, one should remember that the soliton excitations in the 
spiral phase lie in the high-energy part of the spectrum, because the energy of the soliton-likc solutions differs from 
the ferromagnetic one on some finite value (|26|) . Therefore, the soliton energy is higher than the spiral one (|23[) on 
the value proportional to the system size. Nevertheless, these states can be important in strong magnetic field close 
to the saturation value. 

If we expand n z w 1 — \ip\ 2 /2 and take into account the terms ~ \ip\ 2 i(j in Eg . ([29f . then Eq.(|29[) near the border of 
the soliton existence reduces to the modified non-linear Schrodinger equation containing the fourth-order derivative 
term 

M 19V, #V , „ , U + iVK 2, , q „x 

l dT = 4W + (l W + 2iJ ~^^ H * (37) 

Unfortunately, in spite of substantial simplification of the initial equation of motion (|29p . the exact solution of this 
equation is unknown and we have to use numerical calculation. 

We note that Eqs. (j3~lj) , (|33[) are valid for small a <C 1, but for any 7. In particular, for the case A = (7 = — 1) 
they correctly reproduce the known asymptotic for 6 and ip. 

IV. NUMERICAL CALCULATIONS 

We have carried out a numerical analysis of the discrete and the continuum versions of the classical spin model @ 
in the vicinity of the IT point. Both energy functionals Q and (|18|) were minimized numerically over angles 9 n , tp n 
and 9(x),(p(x) with fixed values of magnetization {§1 and momentum ([9]). The periodic boundary conditions for 
soliton and the open boundary conditions for kink excitations were imposed. The numerical calculations showed that 
for small a the difference between the discrete and the continuum models is negligible as expected and it increases 
as a grows. On the investigation of the dynamics of solitons we restrict ourself to small anisotropy a -C 1, when the 
continuum approximation is valid, because in the discrete model with finite a a so-called pinning potential appears 
[2^ , I30I ] and the momentum becomes undefined. So, on default we will present numerical results for the discrete 
classical model (jj), keeping in mind that the continuum approach corresponds to the limit a — > 0. 

We studied the finite size effects for both discrete and continuum models and found that the system size N should 
be taken so that the parameter a^-^N 1. The finite size effects will be discussed in detail in Sec.V. Here we only 
notice that when the relation a^^N ^ 1 is fulfilled the convergence of a solution accelerates exponentially with AT. 

At first we present and discuss the results for small values of a and 7, when the continuum approximation is justified 
and the soliton energy takes the scaling form ([26]) . In order to verify the found scaling equation we plotted the results 
for a fixed k and different values of a, 7, M as E/(s 2 a 3 ^ 4 ) vs. scaling parameters v — a}f A M/s and fi = j/y/a. For 
small a and 7 all data must lie on one curve which actually represents the scaling function f(p, v, k). 

For example, the numerical data for the static solitons (k = 0) and for // = is shown in Fig[3J Here the calculated 
dependencies of the soliton energy on the soliton size M for three values of a — 1CP 2 , 10~ 3 , 10~ 4 are demonstrated 
in axes E/(s 2 a 3 ^ 4 ) vs. v — a 1 ' 4 M/s. As we see all three solid curves lie very close to each other and they rapidly 
converges in the limit a — > 0, so that the curve corresponding to a — 10~ 4 perfectly describes the scaling function 
f(0,v,0). 

The function /(0, v, 0), shown in Fig|3l reaches a maximum at v « 5 and then rapidly saturates, yielding the energy 
of large solitons E ~ 4s 2 a 3 / 4 . For comparison we put also in FigGO similar data for k = and A = 1/4 obtained in 
Ref. [19j] for the quantum spin-1/2 chain of length N — 24 and a — 0.05, which represents a rough estimate for the 
spin- 1/2 scaling function f 1/2(0, v, 0). One can see that the behavior of the scaling functions for the quantum s = 1/2 
and the classical s — > 00 models are very similar, though they have different limits at v — > 00. Thus, we believe 
that the scaling law Eq. (|26[) is valid for quantum model @ with general s and the corresponding scaling functions f s 
behave similar to /1/2 and /. 

It is known [2l[ that for A = there is an equivalence between the energies of the classic static solitons of small size 
and the bound boson complexes of the Bose-Hamiltonian which is mapped from the anisotropic Heisenberg model. 
We expect that this mapping is valid for A ^ as well and it allows to determine the functional form of the energy 
for small solitons near the IT point. As it will be shown in Sec. VI the energy of the M-boson bound state of the 
corresponding Bose-model is 



(38) 
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FIG. 3: Scaled energy E/s 2 o?/ A versus scaled soliton size v = Ma 1 / 4 / s for fi = and k = 0. Solid lines are soliton energies for 
three values of a — 10 -2 , 10 -3 , 10 -4 from top to bottom. Circles (joined by dashed line) are energies of magnon complexes 
of quantum s = 1/2 F-AF chain for N = 24 and a = 0.05 [19(. Dot-dashed line is the energy of boson bound states of model 
(EH). 




—i 1 Ore— I 1 i — 

-2-10 1 2 

FIG. 4: The functions G(/3) for the energy of small static solitons (triangles) and Gb(/3) for the energy of boson complexes 
Eq.JTTJ (solid line). 

where = [ijv 1 ! 3 and the dependence Gb{(3) is shown in FigJU 

We found that the numerical data for f(/j,,v,0) at v <C 1 is perfectly fitted by Eg. ([55)) and the function G(f3) for 
the classical solitons perfectly coincides with Gb{0) (see Fig[J|. Therefore, we believe that the scaling functions for 
solitons and boson complexes are the same at v -C 1 and are given by Eq. ([38|) . However, they are different generally 
(see Fig|3]). In particular, at v 3> 1 the energy of the soliton saturates while the boson complex energy diverges. 

For small solitons with k = according to Eq. (|32|) uj is given by the derivative of the right-hand side of Eq. (|38|) . 
Therefore, the boundary between different regimes of asymptotes given by equation fi 2 = 2—uu (see Eq . (|5rjl) ) transforms 
for small solitons to the relation 3/3 2 = 7G — 2(3G' for j3. Numerical calculation of the function G{0) (FigHJ) indicates 
that there is only one solution (3 = —0.57 of the above relation. This means that the oscillation behavior of the 
solitons of size M exists if /i > — 0.57^ 2 ' 3 and that the first two terms in the expansion in small v (|38[) is not sufficient 
to determine the boundary where solitons disappear. 

The scaling function f(n,v,k) for k = 0, 7r/2, 7t at fj, = 1 is demonstrated in Fig[5] As can be seen in Fig[5] 
the scaling function f(l,v,k) converges at v — > oo to the same finite value for all k, which is natural, because the 
spectrum of large solitons is flat. Another observation followed from FigJS] is that the scaling function for fi = 1 
oscillates substantially stronger than that in the case \i = (FigJ3]). 
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FIG. 5: Scaled energy E/s 2 a 3 ^ versus scaled soliton size v = Ma 1 ' 4 /s for = 1 and for different k. 



As we have stated above the necessary condition for the soliton existence is 



dv 



< U) (fJ,,v) 



(39) 



which means that df /dv < 2 for /i < and df /dv < 2 — /i 2 for /i > for any velocity. 

This inequality imposes certain limitations on the soliton size (or the allowable values of v) with given values fj, and 
k. However, the existence regions in the (v, k) plane can not be obtained directly from that in (v, u>) plane, because 
the analytical relations between v, k and v, u near the IT point are unknown in contrast with the integrable case 
A = 0. Therefore, the stability of the soliton for different values of parameters [i, v and k can be established by 
numerical minimization of the energy functional only. Our calculations show that for /i < the solitons of any size 
and for all k exist and inequality (f39|) is satisfied. However, the situation changes for /i > 0. For example, the solitons 
with = 1 and k = are stable for v > 1 only as it is shown on Figj5] As a matter of fact for fi > only part of a 
whole phase space (a half-strip v > 0, < k < tt) is allowable for solitons. In general, the allowable region in (v, k) 
half-strip has a complicated form and can consist of many disconnected parts. For instance, at fi = 1.5 the static 
solitons with k = are stable in the range 3.5 < v < 8 though for k ^ small solitons exist. 

Though a detailed phase picture can be find only numerically, nevertheless we can make some statements about it. 
Since soliton is a localized object, its energy inevitably saturates at v — > oo to some finite value independent of k. 
This means that for large solitons df /dv — > 0. Therefore, for [i < \/2 all large enough solitons with v > v m i n (/j,) are 
allowable. On the other hand deeply in the spiral phase when /i > ^/2 the allowable region is restricted by v < v max (ii) 
and large solitons do not exist. In fact, we did not find large solitons for /i > \[2 in our calculations. 

It is interesting to study the constraints on the soliton size in the isotropic limit of model © at A > 1/4, when 
a — > and /i — > oo. It can be shown that the function f(/i, v, k) at v -C 1 and any k ^ behaves as (const / v 3 ) , where 
const does not depend on a. In this case condition (|39[) reduces to 7M 2 < const. It means that the maximal possible 
size of soliton grows at 7 — > 0. This fact qualitatively agrees with the observed behavior [TEI, [n| of the stability of 
magnon bound complexes with k = tt in the isotropic quantum s = 1/2 F-AF model at A = 1/4. This is one more 
indication of the resemblance of the classical solitons and the quantum bound complexes in the F-AF model. 

In Figl6] the dependence of the soliton energy on the momentum for fj, = 1 is plotted for three different soliton 
sizes corresponding to v = 3, v = 6 and v = 9. It is interesting to note that similarly to the case A = 0, there is a 
periodical dependence of the energy on the momentum in the continuum model near the IT point. We see that the 
dependence of the energy on the momentum rapidly flattens with the increase of the soliton size. In other words, the 
soliton mass exponentially grows with the soliton size. In this respect the behavior of the soliton spectrum is similar 
to the exactly solvable case A = (see Eq- fTU]) ). However, the dependence E(k) has a very complex form in contrast 
to the case A = 0. As it is shown on Fig[6] /c m i n can be k = 0, k = tt or even intermediate value between them. For 
v = 6 the minimum of E(k) is obtained at k = tt. With the subsequent increase of the soliton size v the form of the 
spectrum alternates: it has the minimum either at k — or at k = tt. However, for v > 10 the spectrum becomes so 
flat that it is difficult to distinguish these cases numerically. The soliton velocity v is zero in both points k = and 

k = TT. 
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FIG. 6: Energy momentum spectra at /i = 1 for v = 3 (dot-dashed line), v = 6 (dashed line) and ^ = 9 (solid line). 




FIG. 7: Shapes of small (z/ = 2) and large (y — 20) solitons for /i = and fc = 0. 




FIG. 8: Shapes of small (y = 2) and large (y = 20) solitons for \i — 1 and fe = 0. 
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FIG. 9: Shapes of small (y = 2) and large (y = 20) solitons for [i — 1 and k — ir. 



We have obtained numerically the soliton shapes for different values of the parameters \x, v and k. The shapes of 
static (k = 0) solitons at fi = for the cases v = 2 and ^ = 20 are demonstrated in Fig(7] The soliton solutions 
oscillate around 9 = (and 9 = ir for large soliton) and these oscillations exponentially decay as moving off the 
domain walls. Numerical calculations show that at approaching to the spiral phase the amplitude of the oscillations 
grows and its damping decreases with the distance from the soliton center. The shape of solitons of two different sizes 
with v = 2 and v = 20 at fi = 1, which is very close to the transition line, is shown in FigjS] Here one can see the 
increase of oscillations in comparison with the case [i = 0. The observed behavior of the soliton shape is in full accord 
with conclusions followed from the analysis of Eq. (f3"o) . 

The shape of solitons changes with the momentum k. This change is more pronounced for small solitons: they 
become narrower and higher. At the same time the shapes of large solitons remain almost the same. As an example 
the shapes of small and large solitons for the momentum k = n and fj, = 1 are demonstrated in Fig|9] As shown in 
FigJ5]the function 6{x) has a cusp at the maximum, which resembles the case A = 0. However, the behavior of the 
function ip(£) drastically differs from the case A = [3l| . In particular, (p(j£) is discontinuous at £ = in contrast 
with that for A = 0, where this function is linear near £ = 0. 

One more illustration of the oscillation behavior of localized excitations in the F-AF model is the kink solution for 
the particular case fi = (A = 1/4). Though we did not find the analytic solution of the non-linear equation (j20|) . we 
found for this particular case some signs that it can admit a solution in a closed form. Namely, we studied numerically 
the kink excitation described by equation (|20[) with 7 = 0. The boundary conditions are 6{— 00) = 0, 9{oo) = 7r and 
6'{—oo) = 9'(oo) = 0. Numerical solution showed with very high accuracy that the asymptotic of the kink solution 
centered at £ = has a form 

6(S -00) = 2sin(l)cos(-2 1/4 £- l)e 2l/4 « 

9{i -> 00) = 7r-2sin(l)cos(2 1/4 £- l)e- 2l/4 « (40) 
The oscillation behavior of the kink asymptotes is in full accord with Eq. ([36]) . 

At iei < 1 

^)-f + 2 1/4 £ (41) 

Besides, the calculated kink energy has a surprisingly simple form 

E = 2s 1 o?l A (42) 

where the numerical factor 2 was found with high precision. Moreover, the contributions to this kink energy from 
different terms in Eq. (flT)|) have simple rational ratios. All these facts together give a hope to obtain the exact solution 
of Eq. (|2"0)) for 7 = in closed form. 

The continuum approximation is valid for a <^ 1. If the parameter a is not very small there are corrections to the 
soliton energy due to a discreteness. The energy of large soliton (^3>l)atA = l/4 obtained numerically is shown in 
Fig llOl As one can see the energy is perfectly described by equation 



E/s 2 » 4a 3 / 4 + 1.12a 5 / 4 



(43) 
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FIG. 10: Dependence of the scaled energy of large solitons on the anisotropy ^fa for A = 1/4. 
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FIG. 11: Functions /(/x) (solid line) and g(fi) (dashed line) describing two leading terms for the energy of large solitons 
(Eq.flHJ). 

which is valid up to a ~ 1. Here the first term is described by the continuum approximation, while the second term 
represents the correction coming from the discreteness of the lattice. The found correction to the soliton energy 
(~ a 5 / 4 ) differs from the correction to the energy of multimagnon complex (~ a) found for the case s = 1/2 in 
Ref.[3. 

We found that similar behavior of the energy of the large soliton takes place for any fi. So, the energy for v ^> 1 
can be written as 

E = s 2 a 3 ^f^) + 5 2 a 5/4 ff (M) (44) 

where f(fJ.) is the saturated values of f((i, v, k) at v — > oo. 

The behaviors of the functions f(fi) and g(n) are shown in FigQT] As one can see /(/z) goes down with the increase 
of /i and becomes negative in the spiral phase at fi > 1.2. At the same time, the correction in Eq. pi)) increases at 
approaching to the spiral phase. 

In the opposed limit of large anisotropy a 3> 1, the asymptotic of the soliton energy is found using a simple 
perturbation theory, which gives 

2 9\ 2 

E/s 2 =4a-8X + 0(a~ 2 ) (45) 

a as 

As follows from this equation the energy is not simply proportional to s 2 in contrast with the case A = (Eq. (|14|) ). 
This means that the dependence of the energy of kink or large soliton is not a universal function of a for general s as 
it is for A = 0. 
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V. FINITE-SIZE EFFECTS. TRANSITION FROM UNIFORM TO LOCALIZED SOLUTION. 



In the preceding sections we investigated the soliton and the kink excitations of model ((4]) in the infinite system. 
Now we study the finite-size effects of these states. An interesting property of the soliton solutions of the classical 
model (01 on finite ring is the existence of the critical value cx.q(L) below which the uniform solution (9(x) — const) 
is realized. At a > ao this state develops into localized soliton. This transition from the uniform to the localized 
solution is similar to the well-known Gross-Pitaevskii transition in the Bose-systems [32| . 

We consider the static soliton on the ring of size L with periodic boundary conditions. We use the continuous 
approach because the finite-size effects are essential at small values of a. In the beginning we study the case A = 0. In 
this case the function 9(x) providing the extremum of the energy functional ([7]) satisfies the Lagrange-Euler equation 



q2q L L 

= asin(20) - aw sin 0, ~~2 <x< ~2 



with boundary conditions 

0'(O) = 9'(±L/2) = (47) 
The Lagrange multiplier lj ensures the condition of a given normalized magnetization 

L/2 

cos0— (48) 

-i/2 

Simple analysis shows that for any given values of n z and L there is a critical value of parameter a = ceo(n z , L) so 
that at a < ao Eq. (j46|) has a single and uniform solution 8(x) = 9q = arccos(n^). The energy of this solution is 

E = Ls 2 a{\ - nl) (49) 



At a > ao another, non-trivial solution of Eq. (|46p appears, which is a precursor of localized soliton solution. In 
particular, for the case n z = the critical value is ao = 2ir 2 /L 2 and the uniform solution at a < ao is 6q = n/2. The 
non-trivial solution at a > a can be expanded in small parameter $ = y/{a — a )/2a <C 1: 

0(0 = tt/2 + 2i? cos £ + 2tf 3 £ sin — cos(3£) (50) 

where £ = (2a) 1 ^ 2 x. 

The energy for this case can be written as 

E 



= I, I < Iq (51) 
E =i -^('-'°) 2 , ^ (52) 



where I = L^fot and lo = n\/2. Therefore, the second derivative of E with respect to I is discontinuous at I = Iq 
(a = a ). 

At / — > oo the solution for the case n z = becomes 

tan (0/2) = exp [V2o~(x + L/4)l (53) 

and E tends to the known results 4s 2 v / 2a. 

A similar behavior of the soliton solution takes place for other values of A. It is straightforward to obtain the value 
ao and the soliton solution near ao for the general case in the same manner as for the case A = 0. For example, for 
A = 1/4 and n z — the critical value is ao = 27r 4 /L 4 and the solution at a <^ ao is 

0(0 = \ + * cos i + L £ sin £ + L co8(30 (54) 



where £ = (8a 



,1/4, 
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FIG. 12: Scaled soliton energy versus finite-size parameter Not 1 / 4 for A = 1/4 and k = 0. Circle is a critical point of soliton 
creation. 



The energy in this case is 



E 



s 2 a 3 / 4 
E 

s 2 a 3 / 4 



= l- 



l, I < l 
3(1 ~l ) 2 



2 1 / 4 >TT 



l>lo 



(55) 



where I = a x ^ 4 L and l = 2 1 / 4 7r . 

The calculated dependence E(l) for A = 1/4 and n z = is shown in FigfT^l As one can see the energy rapidly 
saturates to its asymptotic value at I > Iq. Therefore, in order to reduce finite-size effects in numerical calculations 
near the transition point A = 1/4 one should choose the system size N so that the parameter a^-^N > 10. 

The critical ag for a given magnetization n z can be restored from the case n z = by the relation 



a (n Zl L) 



1-n? 



(56) 



The process of the formation of the kink (domain-wall) occurs in a somewhat different way. The kink is realized in 
open chain with the boundary conditions 9(—L/2) = 0, 0(L/2) — tt. In this case the lowest state at small a is the 
spiral one 



. TT TTX 

ew = - + T 

and its energy remains finite at a = in contrast with the solitons. For example, for A = 1/4 it is 



E 



2 2tt 4 



, aL 



(57) 



(58) 



When a increases the spiral configuration continuously transforms to the kink. In this respect the situation is 
different from the definite and sharp transition from the uniform to soliton-like state studied before. The crossover 
between the spiral and the kink states takes place at I ~ 1. At / — * oo the energy of the kink is a half energy of the 
soliton. 

It is worth noting that the transition of such kind coming from the finite-size effects occurs in the quantum spin 
model ([2]) as well. For very small a, when the corresponding parameter I <C 1, the system is in a uniform state, which 
is the state with a given total projection S z and the maximal total spin S = 5 max = sN. The energy of this state 
is given by Eq.l[49|) with n z = S z /sN as in the classical model. This uniform state transforms continuously to the 
localized state of bound magnon complex when a increases. The crossover between the uniform and the localized 
states takes place at a^^N ~ 1 for A = and at a^-^N ~ 1 for A = 1/4. 
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VI. QUANTUM MODEL. MAPPING TO <5- ATTRACTIVE BOSE MODEL. 



It is known that there is a resemblance between the classical solitons and the quantum magnon complexes at 
A — 21]. In the preceding sections we observed such a resemblance for the classical F-AF model and its quantum 
counterpart with s = 1/2 as well. Here we consider the excitations in the quantum model with general s and compare 
them with the classical solitons. A standard method to treat the quantum spin models is a mapping of the spin- 
Hamiltonian to the Bose-one. Though this method is approximate it gives important indications about the behavior 
of the spin systems. Using the Dyson-Maleev transformation 



S+ = VTsaf { 1 - 



a- cii — s 

afai 



2s 

Sr = V2sa t (59) 
we represent the Hamiltonian ^ in the form 

H = ^2 £kat a k + — v ( k i^ k 2,q)aq +ki a+_ ki a q+k2 a q ^ k2 (60) 

k k±,k2,q 

where a k and a k are conventional bose-operators. 

In the long- wavelength and weakly-coupling limit V{k\, k%,q) = —a and the one-magnon spectrum e k in the vicinity 
of the IT point is 

sfc 4 

£fc = 2sa — s^fk 2 H — — (61) 
In this limit Hamiltonian (|60[) is equivalent to the Bose-model with the attractive S -function interaction 

m / Q2 i QA \ m 

H = 2sma + sJ2 7a + TZ I " 2a E ^-xj) (62) 



j>i=l 



where m denotes number of the Bose particles. For the case A = one can neglect the fourth-order derivative term and 
such model becomes exactly solvable one [331 ] . However, in the vicinity of the IT point the fourth-order derivative term 
in Eq. (|62"]) is important and it destroys the exact solution found in Ref.[33j]. Therefore, we have to use approximations. 

First, we consider the two-boson bound state near the IT point and compare it with the two-magnon bound state 
found in Ref . [l9j . The wave function of two bosons with the total momentum 2k and the energy Ei is ty{x\,x-i) = 
e it(i 2 +ii)^j, 2 _ Xl y The function ijj is determined from the Schrodinger equation 

(2 S7 - 3 S fc 2 )^^ + _ 2a5( x ^( x ) = e^{x) (63) 

where x = — x\ and — 2e k — E2 (Eb > 0) is the binding energy. Here we see that the fourth-order derivative 
term introduces the dependence of the function tp(x) and the binding energy Eb on the total momentum. 

The function ip(x) as well as its first- and second-order derivatives are continuous at x = 0, while the third-order 
derivative satisfies the condition 



s o"ip 
4 'dx 1 



= o^(0) (64) 

-0 

The wave function of the bound state of two particles with 5 -attraction has a form: 

ip{x) = e~ Klxl [t]cos(t]x) + Ksin(T?|x|)] (65) 

where the parameters k and r\ are determined by the above stated conditions for the function tp(x) and its derivatives 
at x = 0. The solution has a form 



« - (S)_ /<« 

V = g) 1/3 vW?M (66) 
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where 





3(/3) 



7-3fc 2 /2 
(2a/s) 2/3 

1 + xA+ 



(4/3/3)= 



1/3 



4/3/3 



1 + v/1 + (4/5/3)3 



1/3 



(67) 
(68) 



The binding energy is 



£6 



In particular, for 7 = (A = 1/4) and k — it is 



2 1 / 3 q 4 / 3 



a 4 / 3 
(2s)V3 



(69) 



(70) 



The wave function (|65[) has one interesting feature. It oscillates when 77 is real and has a form of superposition 
of two decaying exponents otherwise: ip( x ) = e~ Kl \ x \ + ce _K2 ' :E '. The boundary between these regions is defined by 
equation rj — 0, which has a solution (3 = — 1/2 5 / 3 . Thus, the wave function of two coupled bosons oscillates for 
27 > 3k 2 — (a I s) 2 / 3 and does not oscillate for 27 < 3fc 2 — (a/s) 2 ^ 3 . The oscillation region diminishes with the increase 
of the total momentum. Certainly, the oscillatory behavior of the wave function is the effect of the frustration similarly 
to that observed for the solitons in the spin model. 

We can compare Eq. (|70| with the binding energy of two magnons for s = 1/2 found by us in (l9l | 



E b = a 4 / 3 ~ |a 5 / 3 + 0{a 2 ) 



(71) 



It follows from Eqs. (|7T)|) and (|7Tj) that the leading term in Eb at a — » coincides up to numerical factor with the 
two-boson binding energy, but subsequent terms in the expansion of Eb are absent in (|70|) . Similar difference between 
Eb and e takes place for any values of A. For example, for A = it is 



£6 



of 
27 

2s 



4s- 1 
~8?~ 



a 3 + 0(a A ) 



(72) 



As will be shown below, this difference plays a key role at comparison of the multimagnon and multiboson bound 
states. 

The exact wave function of the m-boson state at A = is known for any m [33| . As was shown in Ref. [34| at A = 
the Hartree approximation correctly reproduces the exact energy of the multiboson states (m ^> 1) with zero total 
momentum. Therefore, we expect that this approach gives reliable results for m 3> 1 when A ^ 0. Here we consider 
Bose model (R)2"j) in a parametric regime corresponding to the vicinity of the IT point of model @ . 

The energy functional for zero total momentum in the Hartree approximation is 



E m = 2sma + 771 



dx 



'dcp\ 2 s (d 2 4> 



(73) 



Reseating £ = (ma/s)^ 3 x and — (ma/s) 1 / 6 x(£) transforms the energy functional ([73]) to 



/oo 
-oo 



dx 



1 (d 2 x 



4 V^ 2 



(74) 



where [3 — 7/(ma/s) 2 / 3 . The function 4>(x) satisfies the conventional normalization condition, which after rescaling 
gives the same normalization condition for x(£) 



/oo 
X 2 (C)dC = 1 
-oo 



(75) 
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The Hartree equation comes from the minimization of the energy over x(£). For m 3> 1 it has a form: 



d 2 x ld 4 x „ 3 / 7R x 

P W + *W (76) 



where e is the Lagrange multiplier secured the norma condition (|75|) . 

Unfortunately, the solution of this equation at present is unknown. However, if we assume the existence of a 
localized solution of Eq. (|75|) . then the integral in Eq. JTi)) converges yielding some function of parameter (3. So, the 
energy of m-boson bound complex takes the form 



E rn = 2sma - sm(— ) 4/3 G b {(3) (77) 
s 



The behavior of the function G b (/?) obtained by the numerical minimization of the functional (|74j) is demonstrated 
in FigJU Numerical calculations also showed that there is a critical value of the parameter (3 rj —0.57 separating 
the regions with ((3 > —0.57) and without ((3 < —0.57) oscillations of the function 4>(x). The value (3 — —0.57 is 
the same as was found for the classical spin model, which certificates the distinct correspondence between the boson 
and the classical spin models for small solitons. However, the bound complex of bosons exists in the whole region of 
parameter (3, which means that the instability of solitons studied in Sec. Ill comes from the effects neglecting in the 
boson approach. 

It is interesting to compare Eq. (|77f with the energy of the m-magnon bound state of the quantum model ([2]) with 
s = 1/2 at A = 1/4. The latter has been obtained in Ref. [l9l| and has a form 

E m — ma — Ci m 7/3 a 4/3 + C 2 m 11/3 a 5/3 + ... (78) 

where C\ and Ci are numerical coefficients. 

Similarly to the case m = 2 the leading terms of the expansions of these two energies in small a coincide with 
Eq. (|77[) up to numerical factors. However, energy (|7T[) tends to — oo at m — > oo (collapse phenomenon), while the 
energy of m-magnon bound complex in the spin model are finite at m — ► oo and behaves as E m ~ a 3 / 4 (l9j . Therefore, 
subsequent terms in the energy expansion in a of the spin model are responsible for a short distance repulsion of 
magnons preventing the collapse. 

A comparison of the expression for the energy of the classical soliton of size M (|26|) with a formula for the energy of 
m-boson complex (|77p indicates that the obtained expression for the energy of m-boson complex (|TT[) is a particular 
case of a more general scaling relation (j2l)|) for k = and 

f b ( f i ) v,0) = 2v-v 7 / 3 G b (-^) (79) 

As we discussed in Sec. IV the scaling function of the boson model fb(fi,v,0) at v <C 1 coincides with the scaling 
function of the classical spin model. The comparison of this equation with Eq. ([26l) allows us to assume that Eq. (|79"l) 
represents two first terms of the scaling function of the quantum spin-s model (including the classical limit) in small 
parameter v 

oc 

f s (^,0) = J2^ n/3+1 9n,s(^) (80) 

n=0 

where functions g UtS depend on s. 

Eq. ([75]) reproduces (probably exactly for any s) two first terms in this expansion. But the expansion ([8H)) contains 
infinite number of terms in contrast with Eq. (|T9"|) . It leads to the finite energy of both the classic solitons and the 
quantum bound magnon complexes while the energy of m-boson complex diverges at m — ► oo. 

Summarizing all the above facts, we believe that the bound energy of m -magnon complex for a quantum spin-s 
model ((2]) near the IT point is correctly described by the classical scaling formula (|26|) 



E = s 2 a 3 / 4 f s {n,v,k) (81) 

The scaling functions f s for quantum spin-s case does not coincide with the function / obtained in the classical 
continuum approach, though all of them have very similar behavior and lim f s = f. 
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VII. CONCLUSION 

We studied the soliton excitations in the classical F-AF model with the easy-axis anisotropy. The F-AF model has 
two parameters: the frustration parameter A = | J2I/J1 and the anisotropy a. We found that in a weakly anisotropic 
limit (a <C 1) the behavior of the soliton solutions for small frustration parameter A < 1/4 is qualitatively similar to 
that for the exactly solvable easy-axis XXZ chain (A = case). However, the situation drastically changes near the 
IT point {a = 0, A = 1/4), where the transition from the ferromagnetic to the spiral ground state takes place in the 
isotropic case. In the vicinity of the IT point the corresponding energy functional in the continuum approximation 
qualitatively changes and does not admit the exact solution. The analysis of the derived energy functional allowed 
us to estimate the behavior of the transition line between the ferromagnetic to the spiral ground state in (a, A) plane 
near the IT point. 

We mainly interested in the behavior of the solitons in the vicinity of the IT point. We showed that these localized 
states are separated from the ferromagnetic state by a finite gap. The dependence of the soliton energy (the gap) on 
model parameters near the IT point was established on a base of the scaling arguments. As a result we found that 
the soliton energy is proportional to ~ a 3 / 4 and is expressed by the scaling function (j2l)|) depending on three scaling 
parameters: /1 — (A — 1/4) /y/a, the scaled soliton size v and the momentum k. 

The analysis of the asymptotic solutions of the corresponding equation of motion provided us with the necessary 
conditions of the soliton stability. It was shown that solitons of all sizes and any momentum exist when fi < 0, while 
in the region /1 > there are definite restrictions on the soliton size v and the momentum k. The distribution of 
the allowable values of the soliton parameters in (y, k) plane for /1 > has very complicated form, which can be 
determined numerically only. For example, the static solitons (k = 0) only of the middle size 3.5 < v < 8 are stable at 
fi = 1.5. Nevertheless, some facts about the soliton existence region was ascertained analytically and then confirmed 
by numerical calculations. In particular, small static (k = 0) solitons are unstable for fi > 0, while small solitons with 
non-zero momentum exist for any \x. Large solitons [y ^> 1) exist for \i < y2 and, therefore, they survive in the part 
of the spiral phase 1 < fx < y2. Though the soliton excitations in the spiral phase lie in the high-energy part of the 
spectrum, they can play an essential role in the magnetization processes. 

In the isotropic limit (/x — ► 00) only small solitons with non-zero momentum survive and the maximal allowable size 
of these solitons increases when the frustration parameter tends to the critical value A — > 1/4. Such a dependence of 
the soliton size on the frustration parameter is qualitatively similar to the observation that the size of the multimagnon 
bound complexes with k = -k in the quantum F-AF model grows when A — > 1/4. 

We found that the frustration reveals itself in the oscillating shape of solitons. The amplitude of the oscillations 
grows at approaching to the spiral phase and further inside of it. Generally, at some value of fj, the solution starts to 
oscillate without the decay at infinity, i.e., soliton- like solution with given values of v and k disappears. 

We studied the finite-size effects on the soliton solutions. It was shown that the soliton solution on finite ring 
originates in the uniform (non-localized) state. The transition from the uniform to the soliton-like state occurs at the 
critical value of the anisotropy ao(L) ~ L~ 4 , below which the uniform solution is realized. This finite-size effect is 
similar to the well-known Gross-Pitaevskii transition in the Bose-systems. 

In order to establish a connection between the properties of the solitons and the multi-magnon complexes of the 
quantum counterpart of this model we used the Dyson-Maleev mapping of the spin model to the Bose one. It turned 
out that the dependence of the energy of boson bound complexes on model parameters represents a particular case 
of the found scaling expression for the classical soliton energy. Moreover, the energy of the bound magnon complexes 
for quantum spin- 1/2 model found by us before perfectly coincides with the scaling equation for soliton energy. 
Therefore, we believe that the bound energy of multi-magnon complex for a quantum spin-s model near the IT point 
is characterized by the same critical exponents and by the identical scaling parameters as the soliton energy, though 
the corresponding scaling functions are different for different s. It is known that such a resemblance takes place for 
the Heisenberg ferromagnetic chain with an easy-axis anisotropy. Our study shows that the multi-magnon complexes 
behave substantially as the classical objects for the frustrated model as well. 

For future it would be interesting to study the behavior of the classical spin model deeply in the spiral phase. We 
believe that it can help in understanding of unusual magnetization processes and shed light on the behavior of the 
multi-magnon complexes at high magnetic fields. 
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